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BACKGROUND 



1. Field of the Invention 

This invention relates to adaptive filtering. In particular, the invention relates to 
adaptive filtering with tap leakage. 



2. Description of Related Art 

Adaptive filters may be used in a number of applications such as adaptive 
equalization in digital communications systems. Adaptive filters are typically 
implemented as a non-recursive filter having N taps. For example, the adaptive filter may 
10 be a fractionally spaced equalizer having tap weights spaced a fraction of a symbol. The 
fractional spaced equalizer provides better performance than conventional symbol 
synchronous equalizers in the presence of severe linear distortion. The adaptive filter 
may have N equalizer coefficients that are updated or adjusted according to the input 
signal. The adaptive filter prevents the buildup of large coefficient values by 
15 systematically adjusting the magnitudes of all the equalizer tap weights or equalizer 

coefficients. However the systematic adjustment may cause output overflow, leading to 
substantially degraded performance. 

A technique to solve this problem includes preventing the buildup of large 
coefficient values through systematically "leaking" or decreasing the magnitudes of the 
20 equalizer tap weights. This technique has a number of disadvantages. First, it requires an 
additional multiplication for each tap-update. The number of additional operations 
becomes significant for large number of taps. Second, adding tap-leakage computations 
to pre-existing hardware blocks may be difficult, considering the number and complexity 
of additional required operations. 

25 Therefore, there is a need to have an efficient technique to reduce the 

computations of the adaptive filter with tap leakage. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The features and advantages of the present invention will become apparent from 
the following detailed description of the present invention in which: 

Figure 1 is a diagram illustrating a system in which one embodiment of the 
invention can be practiced. 

Figure 2 is a diagram illustrating an adaptive filter shown in Figure 1 according to 
one embodiment of the invention. 

Figure 3 is a diagram illustrating a tap-leakage generator shown in Figure 1 
according to one embodiment of the invention. 

Figure 4 is a diagram illustrating an error filter shown in Figure 3 according to one 
embodiment of the invention. 

Figure 5 is a flowchart illustrating a process for adaptive filtering with tap leakage 
using error filtering according to one embodiment of the invention. 
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DESCRIPTION 

The present invention is a technique to implement a tap-leakage computation for 
adaptive filters. The complexity of the tap-leakage computation is reduced from N 
operations, where N is the number of taps in the adaptive filter, to a computation which is 
independent of filter length. In one embodiment of the invention, a tap-leakage generator 
has an error filter and an updater. The error filter is coupled to an adaptive filter having N 
taps to filter a decision error provided by the adaptive filter using a leakage factor. The 
updater updates N equalizer coefficients to the N taps using the filtered decision error. 
The updater receives the N equalizer data from the N taps. 

In the following description, for purposes of explanation, numerous details are set 
forth in order to provide a thorough understanding of the present invention. However, it 
will be apparent to one skilled in the art that these specific details are not required in 
order to practice the present invention. In other instances, well-known electrical 
structures and circuits are shown in block diagram form in order not to obscure the 
present invention. 

The present invention may be implemented by hardware, software, firmware, 
microcode, or any combination thereof. When implemented in software, firmware, or 
microcode, the elements of the present invention are the program code or code segments 
to perform the necessary tasks. A code segment may represent a procedure, a function, a 
subprogram, a program, a routine, a subroutine, a module, a software package, a class, or 
any combination of instructions, data structures, or program statements. A code segment 
may be coupled to another code segment or a hardware circuit by passing and/ or 
receiving information, data, arguments, parameters, or memory contents. Information, 
arguments, parameters, data, etc. maybe passed, forwarded, or transmitted via any 
suitable means including memory sharing, message passing, token passing, network 
transmission, etc. The program or code segments may be stored in a processor readable 
medium or transmitted by a computer data signal embodied in a carrier wave, or a signal 
modulated by a carrier, over a transmission medium. The "processor readable medium" 
may include any medium that can store or transfer information. Examples of the 
processor readable medium include an electronic circuit, a semiconductor memory 
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device, a ROM, a flash memory, an erasable ROM (EROM), a floppy diskette, a compact 
disk (CD-ROM), an optical disk, a hard disk, a fiber optic medium, a radio frequency 
(RF) link, etc. The computer data signal may include any signal that can propagate over a 
transmission medium such as electronic network channels, optical fibers, air, 
5 electromagnetic, RF links, etc. The code segments maybe downloaded via computer 
networks such as the Internet, Intranet, etc. 

It is noted that the invention may be described as a process which is usually 
depicted as a flowchart, a flow diagram, a structure diagram, or a block diagram. 
Although a flowchart may describe the operations as a sequential process, many of the 
10 operations can be performed in parallel or concurrently. In addition, the order of the 
operations may be re-arranged. A process is terminated when its operations are 
completed. A process may correspond to a method, a function, a procedure, a subroutine, 
a subprogram, etc. When a process corresponds to a function, its termination corresponds 
to a return of the function to the calling function or the main function. 

15 Figure 1 is a diagram illustrating a system 100 in which one embodiment of the 

invention can be practiced. The system 100 includes an adaptive filter 110 and a tap- 
leakage generator 120. 

The adaptive filter 110 receives an input signal and generates an output sequence. 
The adaptive filter 110 may be used in a number of applications such as adaptive 

20 equalization in a digital communications systems. The adaptive filter 1 10 is typically 

implemented as a non-recursive filter having N taps. An adaptive equalizer can be either 
symbol-spaced or fractionally spaced. A fractionally spaced equalizer has tap weights 
spaced a fraction of a symbol internal apart. The fractionally spaced equalizer provides 
better performance than the conventional symbol-spaced equalizer in the presence of 

25 severe linear distortion. The adaptive filter 1 10 is described here as a symbol-spaced 

equalizer, but the invention is valid for both symbol and fractionally spaced filters. The 
adaptive filter 110 has N equalizer coefficients that are updated or adjusted according to 
the input signal. The adaptive filter 110 prevents the buildup of large coefficient values 
by systematically adjusting the magnitudes of all the equalizer tap weights or equalizer 
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coefficients. The adaptive filter 110 generates N equalizer data and a decision error to the 
tap-leakage generator 120. 

The tap-leakage generator 120 updates the N equalizer coefficients used by the 
adaptive filter 110 using the N equalizer data and the decision error provided by the 
5 adaptive filter 110. The tap-leakage generator 120 uses an efficient procedure that 
employs less computations than the traditional update algorithm with tap-leakage. 

Figure 2 is a diagram illustrating the adaptive filter 110 shown in Figure 1 
according to one embodiment of the invention. The adaptive filter 110 includes a 
sampler 210, N delay elements 220 0 to 220n-i, N multipliers 230 0 to 230 N .i, an adder 240, 
10 a thresholder 250, and an adder 260. 

The sampler 210 samples the input signal r(t) into a discrete-time sequence of 
input data. The sampling frequency is in accordance to the frequency of the input signal 
meeting the sampling Nyquist criteria. In one embodiment, this sampling frequency is at 
an integer multiple of the symbol rate corresponding to modem applications. The delay 

15 elements 220 0 to 220n-i each delays the input sequence by one sample time period. The 
sample time period is the sampling time interval of the sampler 210. The delay elements 
220o to 220n-i may be implemented by a shift register, a number of flip-flops connected in 
cascade, or any other suitable elements. When implemented by software, the delay 
elements 220o to 220n-i may be implemented an array of elements. Each of the delay 

20 elements generates an equalizer data R n . For N delay elements, there are N equalizer data 
R n ° to R n N_1 . 

Each of the multipliers 230 0 to 230n-i multiplies one of the equalizer coefficients 
D n ° to D n N_1 provided by the tap-leakage generator 120 with the corresponding equalizer 
data. For example, the multiplier 230 0 multiplies the equalizer coefficient D„° by the 
25 equalizer data R n °, D n k with the equalizer data R n k , etc. The N equalizer coefficients D n ° 
to D n N_1 are provided by the tap-leakage generator 120. 

The adder 240 adds all the products generated by the multipliers to provide a sum 
S n . The thresholder 250 thresholds the sum S n to provide the output sequence A n . The 
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adder 260 adds the sum S n and the output sequence to provide a decision error E n at each 
sampling time. 

Figure 3 is a diagram illustrating the tap-leakage generator 120 shown in Figure 1 
according to one embodiment of the invention. The tap-leakage generator 120 includes 
5 an error filter 310 and an updater 320. 

The error filter 310 filters the decision error E n from the adaptive filter 110 using 
a leakage factor a. The error filter 310 generates a filtered decision error h n at each 
sampling time. 

The updater 320 updates the vector of N equalizer coefficients D n ° to D n N_1 to the 
10 N taps using the filtered decision error h n . The updater receives the N equalizer data from 
the N taps of the adaptive filter 110. The updater 320 includes N coefficient updaters 
3 3 0 0 to 330n_i. Each of the N coefficient updaters 3 30 0 to 330^ generates a 
corresponding updated equalizer coefficient. For example, the coefficient updater 330k 
generates the equalizer coefficient D n k . Each of the N coefficient updaters 330o to 330n-i 
15 uses a corresponding equalizer data. All of the N coefficient updaters 330o to 330m.i use a 
same adaptive step size parameter u. and the filtered decision error h n . 

In the conventional tap-leakage procedure, the updated equalizer coefficients are 
determined by the following equation: 

C n+ i = a*C n -n*E n *R n (1) 

20 where: 

C„ is the vector of equalizer coefficients at time n, Cn T = [C n ° Cn 1 . . . C™~ 1 ] 
where T denote vector transposition. 

R„ is the vector of equalizer data at time n, Rn T = [R n ° R n J . . . R n N_1 ] 
E n is the decision error at time n 
25 (X is the adaptive step size parameter 

a is the leakage factor, 0.0 < a < 1.0. 
* is the multiplication operation. 
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The conventional tap leakage procedure includes the multiplication of a in each of 
the update path for the equalizer coefficients. The z-transform of the update path can be 
written as: 

C^z^a^-^E^R, (2) 

5 or: 

C n = -|u*E n *R n z- 1 /(l-a*z 1 ) (3) 

The updater 320 removes the multiplication of a from each of the update paths. 
In addition, the error filter 3 10 is used to incorporate the leakage factor a. For the tap- 
leakage generator 120, the equalizer coefficients are now defined as where D„ T = [D n ° 
10 Dn 1 . . . D n N_1 ]. The least mean square (LMS) tap update for the leakage generator 120 is 
written as: 

Dn+rB.-A*!. ( 4 ) 
where h n is the filtered decision error. 
The Z- transform of the filtered decision error, h n is: 
15 h n = E n *h(z) (5) 

Taking the z-transform of equation (4) to obtain: 

D n *z 1 = oc*D n - u*E„*H(z)*R n (6) 

or: 



_ -u,*E n R„z / s 
D n =— " 2— £ * H(z) 



1-z- 1 



(7) 



20 Setting the coefficients C n equal to D n and solving for H(z) to obtain: 

H(z) = ( 1-2 1 ) (8) 
l-a*z 

The block H(z) is the error filter 310 which is a first order high pass filter (HPF) 
with a zero at zero frequency or direct current (DC) frequency and a pole determined by 
the leakage factor a. 

25 For C n to be identically equal to D n , the transfer function from E n to C n (and also 

E n to D n ) needs to be linear and time- variant. This is not truly the case because R n is not 
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limited to a constant value. However, since the LMS algorithm stochastically estimates 
the gradient of the data-error product, only the time-average of the data-error product is of 
concern, and being strictly linear, time-invariant is not required. 

The updater 320 therefore eliminates N multiplications due to a and the 
5 computation for the error filter 3 1 0 is reduced to only three computations. This 

represents a significant saving when N is large as is typical in adaptive equalization. 

Each of the coefficient updater 330k includes a first multiplier 342k, a second 
multiplier 344k, an adder 346k and a delay element 348k. The first multiplier 342k 
multiplies a corresponding one of the equalizer data R n k with the filtered decision error h n 

10 to provide a first product. The second multiplier 344k multiplies the first product with 
the adaptive step size jo, to provide a second product. The delay element 348k delays the 
updated equalizer coefficient D n k to provide a delayed coefficient. The subtracter 346k 
subtracts the second product from the delayed coefficient to provide the updated filtered 
coefficient D„ k . As seen, the multiplication of the leakage factor a is eliminated. The 

15 computation of the updated equalizer coefficient for tap j is given by the following 
equation: 

T>J = Dj-v.*h n *R a { (9) 

Figure 4 is a diagram illustrating the error filter 310 shown in Figure 3 according 
to one embodiment of the invention. The error filter 310 is essentially a first-order high 
20 pass filter and can be implemented by a number of methods. In one embodiment, the 
error filter 310 includes first and second computing elements 410 and 420. 

The error filter 310 filters the decision error E„ to generate the filtered decision 
error h n . The filtered decision error h n is used by all the coefficient updaters 330 0 to 330n. 
! as shown in Figure 3. The difference equation for the error filter 310 can be obtained by 
25 taking the inverse Z- transform of equation (9): 

h n = E n -E n .i + a*h n _i (10) 

where E n _i and h n _i are the delayed versions of the decision error E n and the filtered 
decision error h n , respectively. 
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The first computing element 410 computes the difference d n = E n - E„_i and the 
second computing element 420 computes h n = d n + a*h n .i. 

The first computing element 410 includes a delay element 412 and a sub tractor 
414. The delay element 412 delays the decision error E n by one sample to generate the 
5 delayed decision error En_i. The subtracter 414 subtracts the delayed decision error E n _i 
from the decision error E n to provide the error difference d n . 

The second computing element 420 includes an adder 422, a delay element 424, 
and a multiplier 426. The delay element 424 delays the filtered decision error h n to 
provide a delayed output h n _i. The multiplier 426 multiplies the leakage factor with the 
10 delayed output to generate a product P = a*h n .i. The adder 422 adds the error difference 
d n to the product P to generate the filtered decision error h n . 

Figure 5 is a flowchart illustrating a process 500 for adaptive filtering with tap 
leakage using error filtering according to one embodiment of the invention. It is noted 
that the START may be located anywhere in the process 500. 

15 Upon START, the process 500 generates equalizer data R„ using the N tap delay 

elements (Block 510). For N taps, the vector of equalizer data R„ includes R n ° to R n " . 
The generation of the equalizer data R n may be performed by shifting the elements one 

sample while acquiring the new sample. Next, the process 500 computes the sum S n = X 

R n ' * Dj, where i = 0, . . ., N-l (Block 520). Then, the process 500 thresholds the sum S„ 
20 to generate the output sequence A n (Block 530). Next, the process 500 computes the 

decision error E n = A„ + S n (Block 540). Then, the process 500 filters the decision error 
E„ by computing the filtered decision error h n = E n - E n .! + a*h n _i in accordance to 
equation (10) (Block 550). The filtered decision error h n will be used to update the 
equalizer coefficients D n . 

25 Next, the process 500 updates the equalizer coefficients p_n using the filtered 

decision error h n (Block 560). The updated equalizer coefficients D n are computed 
according to equation (9). 
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Next, the process 500 determines if the termination condition has been met (Block 
570). This termination condition may be based on a maximum number of iterations, a 
maximum time period, or when the system stops working. If so, the process 500 is 
terminated. Otherwise, the process 500 returns to Block 510. 

While this invention has been described with reference to illustrative 
embodiments, this description is not intended to be construed in a limiting sense. Various 
modifications of the illustrative embodiments, as well as other embodiments of the 
invention, which are apparent to persons skilled in the art to which the invention pertains 
are deemed to lie within the spirit and scope of the invention. 
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